Dissipative hydrodynamics in 2+1 dimension 
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In 2+1 dimension, we have simulated the hydrodynamic evolution of QGP fluid with dissipation 
due to shear viscosity. Comparison of evolution of ideal and viscous fluid, both initialised under the 
same conditions e.g. same equilibration time, energy density and velocity profile, reveal that the 
dissipative fluid evolves slowly, cooling at a slower rate. Cooling get still slower for higher viscosity. 
The fluid velocities on the otherhand evolve faster in a dissipative fluid than in an ideal fluid. The 
transverse expansion is also enhanced in dissipative evolution. For the same decoupling temper- 
ature, freeze-out surface for a dissipative fluid is more extended than an ideal fluid. Dissipation 
produces entropy as a result of which particle production is increased. Particle production is in- 
creased due to (i) extension of the freeze-out surface and (ii) change of the equilibrium distribution 
function to a non-equilibrium one, the last effect being prominent at large transverse momentum. 
Compared to ideal fluid, transverse momentum distribution of pion production is considerably en- 
hanced. Enhancement is more at high pt than at low pr- Pion production also increases with 
viscosity, larger the viscosity, more is the pion production. Dissipation also modifies the elliptic 
flow. Elliptic flow is reduced in viscous dynamics. Also, contrary to ideal dynamics where elliptic 
flow continues to increase with transverse momentum, in viscous dynamics, elliptic flow tends to 
saturate at large transverse momentum. The analysis suggest that initial conditions of the hot, 
dense matter produced in Au-|-Au collisions at RHIC, as extracted from ideal fluid analysis can be 
changed significantly if the QGP fiuid is viscous. 

PACS numbers: 47.75.+i, 25.75.-q, 25.75.Ld 



I. INTRODUCTION 

Lattice QCD predicts that under certain conditions 
(sufficiently high energy density and temperature), or- 
dinary hadronic matter (where quarks and gluons are 
confined), can undergo a phase transition to a decon- 
fined matter, commonly known as Quark Gluon Plasma 
(QGP). Nuclear physicists are trying to produce and 
detect this new phase of matter at RHIC, BNL. Re- 
cent Au+Au colhsions at RHIC, indicate that dense, 
color opaque medium of deconfined matter is created in 
very central collisions 1, 2, 3, |j|- It is now also un- 
derstood that the quarks and gluons strongly interact, 
giving rise to the notion of sQGP (strongly interacting 
QGP). The experimental data are successfully analysed 
in a ideal fluid dynamic model p\. Hydrodynamic evo- 
lution of ideal QGP, thermalised at Ti— 0.6 fm, with cen- 
tral entropy density of 110 fm~^, or energy density of 
35 GeV/ fm? can explain a large volume of RHIC data, 
the Pt spectra of identified particles, the elliptic flow etc 
[3. However, experimental data do show deviation from 
ideal behavior. The ideal fluid description works well in 
almost central Au-|-Au collisions near mid-rapidity at top 
RHIC energy, but gradually breaks down in more periph- 
eral collisions, at forward rapidity, or at lower collision 
energies Q, indicating the onset of dissipative effects. 
To describe such deviations from ideal fluid dynamics 
quantitatively, requires the numerical implementation of 
dissipative relativistic fluid dynamics. 
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Eckart [3 and Landau and Lifshitz [g formulated the 
theory of dissipative relativistic fluid. Their theories are 
called first order theories and suffer from the problem of 
causality; the signal can travel faster than light. 1st or- 
der theories assume that the entropy 4-current contains 
terms upto linear order in dissipative quantities. This 
restriction results in parabolic equations, which leads to 
causality problem. The causality problem is removed in 
the 2nd-order theories, formulated by Israel and Stewart 
[3. In 2nd order theories, dissipative fluxes are treated 
as (extended) thermodynamic variables and the linear re- 
lation between entropy 4-current and dissipative quanti- 
ties is extended to include quadratic terms. The entropy 
4-current contains terms upto 2nd order in dissipative 
forces. The resulting equations are hyperbolic in nature 
and the causality problem is removed. Naturally 2nd or- 
der theories are more complicated. In addition to the 
usual energy-momentum conservations equations, relax- 
ation equations for dissipative fluxes are required to be 
solved simultaneously. 

Though the theories of dissipative hydrodynamics 
0, la, I3 has been known for more than 30 years , signif- 
icant progress towards its numerical implementation has 
only been made very recently [10|, [111, [ij, [Tj, [Tjl . Ear- 
lier attempts were restricted to simple one-dimensional 
Bjorken expansion [13 ■ Since the hot dense matter, pro- 
duced in RHIC collisions, undergoes large transverse ex- 
pansion, such models are not of much help in extracting 
initial conditions of the fluid from the experiment. Re- 
cently Teaney 10] solved the hydrodynamic equations in 
a 1st order dissipative theory He used the blast wave 
model and calculated the corrections to thermal distri- 
bution functions considering shear viscosity only. Blast 



wave models lack dynamics. Freeze-out surface is pa- 
rameterised and one fit the freeze-out parameters with 
experimental observables. Information about the initial 
conditions of the hot dense matter is not available in 
blast wave models. In 12i,2M, 2nd order theories for dis- 
sipative fluid, is solved for QGP fluid. While the mod- 
els are dynamic, (energy-momentum conservation equa- 
tions are solved) energy-momentum conservation equa- 
tions are solved in 1+1 dimension, assuming longitudinal 
boost-invariance and cylindrical symmetry. Thus these 
models do not give any information on the most sensi- 
tive experimental observable, the elliptic flow. Recently 
in |l4| explicit equations describing the space-time evo- 
lution of non-ideal relativistic fluid, undergoing boost- 
invariant longitudinal and arbitrary transverse expan- 
sion, are given. Equations are written for both the 1st 
order and 2nd order theory. 

At the Variable Energy Cyclotron Centre, Kolkata, 
we have developed a numerical code (AZHYDRO- 
KOLKATA) to solve, both the 1st order and 2nd or- 
der dissipative hydrodynamics in 2+1 dimension. In the 
present paper we will present AZHYDRO-KOLKATA re- 
sults for the 1st order dissipative hydrodynamics. Results 
for 2nd order dissipative hydrodynamics will be presented 
in a later publication. As mentioned earlier, 1st order 
theories are acausal, signal can travel faster than light. 
In one dimension, 1st order theories can be solved analyt- 
ically. Analytical solutions indicate that at early time, an 
unphysical reheating of the fluid can occur +Z, +o| . Even 
though in the present simulation, we donot find any indi- 
cation of early reheating, nevertheless, we maintain that 
only 2nd order results, which do not have the unphysical 
causality the problem are more reliable. In this paper we 
will be concerned mainly with the effect of viscosity on 
fluid evolution and its affects the particle production; px 
distribution and elliptic flow. No attempt will be made 
to explain experimental data. 

Paper is organised as follows: in section II, we describe 
briefly the hydrodynamic equations needed for 1st order 
dissipative hydro-dynamics. In section III, we describe 
the equation of state, shear viscosity coefhcient and ini- 
tial conditions used in the present study. With dissipa- 
tion, equilibrium distribution function is changed. Cor- 
rection to equilibrium distribution function due to non- 
equilibrium effects are described in section IV. We have 
studied particle (pion) production in the model. The 
relevant equations for pion production with the equilib- 
rium distribution function and its correction due to non- 
equilibrium effects, are discussed in section V. Results of 
the numerical simulations are shown in section VI. Lastly, 
in section VII, summary and conclusions are drawn. 



II. 1ST ORDER DISSIPATIVE FLUID 
DYNAMICS 

Relativistic dissipative hydrodynamics and associated 
equations in 2+1 dimensions has been discussed in de- 



tail in ref . [lj| . Any fluid dynamical approach starts from 
the conservation laws for the conserved charges and for 
energy-momentum. For a singly charged fluid, the con- 
servation laws are: 



d^N" = 



(2.1) 
(2.2) 



It must also ensure the second law of thermodynamics 



d^S^ > 0, 



(2.3) 



where S"^ is the entropy current. 

In the present paper we restrict ourselves to central ra- 
pidity region, where the QGP fluid is essentially baryon 
free. We thus neglect Eg 12. II To keep the calcula- 
tions simple, we consider the most important dissipative 
term, the shear viscosity and neglect the other dissipa- 
tive terms, e.g. heat conduction, bulk viscosity. For a 
baryon free fluid, the heat conduction is zero. Bulk vis- 
cosity is zero for QGP (point particles). Bulk viscosity 
will be non-zero in the hadronic phase but presently we 
neglect it. 

With the help of hydrodynamic 4-velocity u^ (nor- 
malised as u^Ufj, = 1) and projector A^'' = 5^'' — u^^u"^, 
with only shear viscosity as the dissipative flux, in Lan- 
dau's frame, the energy momentum tensor, and the en- 
tropy 4-current, can be decomposed as. 
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where e = u^jT^^^u^ is the energy density, p is the lo- 
cal pressure, tt'^" is the non-ideal part of the energy- 
momentum tensor, the stress tensor due to shear vis- 
cosity, s = u^S'^ is the entropy density and $'^ is the 
entropy flux due to dissipation. 

In the first order theories, the shear stress tensors are 
written as. 



TTf"" = 277 V' 



(2.6) 



where -q is the shear viscosity coefficient and V^^u"^^ is 
a trace-less symmetric tensor defined as. 



V<^u''> 



- (A^'^A'^^+A'^'A''^) - ^A^^^A" 



(2.7) 

Viscous pressure tt'''^ is symmetric (tt'''^ = tt"^^), trace- 
less {'k'^ = 0) and transverse to hydrodynamic velocity, 
(m^tt^" = 0). The 16-component tt^'' has only 5 inde- 
pendent components. As mentioned in the beginning, we 
have solved the equations with the assumption of longitu- 
dinal boost-invariance. With boost-invariance the num- 
ber of independent shear-stress tensors further reduces 
to 3. 



Heavy ion collisions are best described in terms of 
proper time r = Vi^ — z'^ and rapidity 77^ = 5 In |^ 
(we use the subscript s to distinguish spatial rapidity 
from the viscous coefficient r\). In (r, x, j/, r\s) coordinates, 
with longitudinal boost-invariance, the hydrodynamic 4- 
velocity can be written as, 



(u^u^'a^'a''=) 



(7j_,7_Lt;j,,7_Li;y,0), 



(2., 



with 7_L — ^1 \ h- 



I — v^. Explicit equations for 

energy-momentum conservation, in (r, x, y. ris ) coordi- 
nate system has been developed in ref.|lj|- Here we 
rewrite the results in a form suitable for numerical al- 
gorithm. The energy-momentum conservation equations 
are. 



dr (f^^) + d, (f --iJ,) + dy (f ^X) =-ip + r^TT""), (2.9) 

a,(f^-) +5,(f-"«.) +9,(f-x) ^-d,{p + n---n--v,)-dy{n-y-n-^vy), (2.10) 

dr(T^y) + a,(T^^i>,) + dy(f^yvy') = -d,{n-y - n^yy,) - dyip + ^yy - n-yvy), (2.11) 

where Vx = T'^^/T'^^ and Vy — T'^y jT'^'^ ^ and we have used the notation "tilde" to represent quantities multiplied by 
the factor t, p = rp and similarly T^^ = tT^^ . We note that unlike in ideal fluid, in viscous fluid dynamics conservation 
equations contain additional pressure gradients containing the dissipative fluxes. Both T"^^ and T'^y components of 
energy-momentum tensors now evolve under the influence of additional pressure gradients. 

In 1st order theory, the shear stress tensor components required in the preceding equations are. 
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'r" = 2r7[^(7i-l) + a.7±-^^(7i)], 

TT-- = 24~-dxh^Vx) ~ ^Di^l) + ^(1 + 7i«D], 
Tryy = 2rj[~dxh±vy) - ^D(li) + ^(1 +7i«^)] 



(2.12) 
(2.13) 
(2.14) 
(2.15) 

(2.16) 
(2.17) 



r 



where D — u^d^ is the convective time derivative, 

D = -f^{dr+uA+Vydy), (2.18) 

and 9 is the local expansion rate, given by, 

^ = — + drj± + dx{vxj±) + dy{vyj^) (2.19) 

T 

Given an equation of state, if energy density (e) and 
fluid velocity {v^ and Vy) distributions, at any time r, are 
known, Eas l2.9l2.1Ul and l2.11l can be integrated to obtain 
e,Vx and Vy at the next time step r^+i. While for ideal 
hydrodynamics, this procedure works perfectly, viscous 
hydrodynamics poses a problem that shear stress-tensor 
components contains time derivatives, 9r7J_, drU^, drU^ 



etc. Thus at time step r^ one needs the still unknown 
time derivatives. In 1st order theories, this problem is 
circumvented by calculating the time derivatives from the 
ideal equation of motion. 



Du^" = 



Wp 



De = -(£ + p)V^w^. 



(2.20) 
(2.21) 



With the help of these two equations all the time 
derivatives can be expressed entirely in terms of spatial 
gradients. 1st order theories are restricted to contain 
terms at most linear in dissipative quantities. Neglect of 
viscous terms can contribute only in 2nd order correc- 
tions, which is neglected in 1st order theories. 



III. EQUATION OF STATE, VISCOSITY 
COEFFICIENT AND INITIAL CONDITIONS 

A. Equation of state 

One of the most important inputs of a hydrodynamic 
model is the equation of state. Through this input, the 
macroscopic hydrodynamic models make contact with 
the microscopic world. In the present calculation we 
have used the equation of state, EOS-Q, developed in 
ref.[3- It is a two-phase equation of state. The hadronic 
phase of EOS-Q is modeled as a non-interacting gas of 
hadronic resonance. As the temperature is increased, 
larger and larger fraction of available energy goes into 
production of heavier and heavier resonances. This re- 
sults into a soft equation of state, with small speed of 
sound, c^ w 0.15. With increasing temperature, the 
available volume is filled up with resonances and the 
hadronic states starts to overlap, and microscopic de- 
grees of freedom are changed from hadron to deconfined 
quarks and gluons. The QGP phase is modeled as that of 
a non-interacting quark (u,d and s) and gluons, confined 
by a bag pressure B. Corresponding equation of state, 
p = \e — |B is stiff with a speed of sound c^ = i. The 
two phases are matched by Maxwell construction at the 
critical temperature, Tc = IQAMeV, adjusting the Bag 
pressure B^/^=230 MeV. As discussed in 0, ideal hydro- 
dynamics explain a large volume of RHIC Au-|-Au data 
with EOS-Q. 



B. Shear viscosity coefficient 



which is 2 times larger than the perturbative estimate. 
In the present paper, we treat the shear viscosity as a 
parameter of the model. To demonstrate the effect of 
viscosity on flow and subsequent particle production, we 
use both the perturbative and ADS/CFT estimate of vis- 
cosity. 

1st order theories are acuasal. As mentioned ear- 
lier, unphysical effects like reheating of the fluid, early 
in the evolution, can occur. In one dimension energy- 
momentum conservation equation can be solved analyti- 
cally. If initial fluid temperature is T^ at initial time r^, 
for constant ri/s, the fluid temperature at time r can be 
obtained as 16], 



T{t) = T, 
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inTi s 

For early times, r < r^nax, 



n 






(7) 
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the solution shows an unphysical reheating. The unphysi- 
cal reheating is minimised if Tmax is small or rj/ s « TiTi. 
As will be explained below, we have used an initial time 
Ti=0.6 fm/c and initial temperature of the fluid, Ti=0.35 
GeV . For both the values of viscosity, 77/s << TiTi, the 
unphysical reheating is minimised. 

Shear viscosity can also be expressed in terms of sound 
attenuation length, Fg, defined as. 



Shear viscosity coefficient (77) of dense nuclear (QGP or 
resonance hadron gas) is quite uncertain. In perturbative 
regime, shear viscosity of a QGP is estimated jl7ikiil8l]. 



77 = 86.473 



1 T3 



(3.1) 



With entropy of QGP, s = 37f^T-^ and a^ «0.5, the 
ratio of viscosity over the entropy, in the perturbative 
regime is estimated as, 



F. 



4t) 
3 

sT 



(3.6) 



Tg is equivalent to mean free path and for a valid hy- 
drodynamic description F^/r << 1, i.e. mean free path 
is much less than the system size. With the present 
choice of equilibration time and temperature, both for 
ADS/CFT and perurbative estimate of viscosity, at ini- 
tial time, Tg/r is much less than unity and hydrodynam- 
ics remains a valid description. At later time the validity 
condition becomes even better. 



'-) w 0.135, 

^ S / pert 



(3.2) 



However, QGP produced in nuclear collisions is non- 
perturbative. It is strongly interacting QGP. Recently, 
using the ADS/CFT correspondence [l9l l2(]||. shear vis- 
cosity of a strongly coupled gauze theory, N=4 SUSY 
YM, has been evaluated, rj = jN^T^ and the entropy is 

given by s = ^^N^T^. Thus in the strongly coupled field 
theory. 



(- 



= -^ « 0.08, 

S / ADS/CFT 47r 



(3.3) 



C. Initial conditions 

As discussed earlier, ideal hydrodynamics has been 
very successful in explain a large volume of data in 
RHIC 200AGeV Au-fAu colhsions "F]. In the present 
demonstrative calculations, we have used the similar ini- 
tial conditions as in ref. 5). Details of the initial con- 
ditions can be found in [5|]. We just mention that in 
ref. 5] , initial transverse energy is parameterised geomet- 
rically. At an impact parameter b, transverse distribution 
of wounded nucleons Nwn{x, y, b) and of binary NN col- 
lisions Nbc{x, y, b) to are calculated in a Glauber model. 



A collision at impact parameter b is assumed to contain 
25% hard scattering (proportional to number of binary 
collisions) and 75% soft scattering (proportional to num- 
ber of wounded nucleons). Transverse energy density 
profile at impact parameter b is then obtained as, 



e{x, y, b) = eo(0.75 x Nwn{x, y, b) + 0.25 x Nbc{x, y, &)) 

(3.7) 

The parameter Sq and the initial equilibration time r^ 
are fixed to reproduce the experimental transverse mo- 
mentum distribution of pions in central Au+Au colli- 
sions. STAR and PHENIX data are fitted to obtain 
initial equilibrium time Ti=0.6 fm and central entropy 
density of s = 110/m~^. This corresponds energy den- 
sity of the fluid as 25 GeV/fm^, or initial temperature 
of 350 MeV. Apart from initial energy density, initial 
velocity distribution is also required in hydrodynamic 
calculation. In the present calculation it is assumed 
that the at the initial time Ti, fluid velocities are zero, 
Vx{x,y) =Vy{x,y) = 0. 

In dissipative hydrodynamics, additionally, initial con- 
ditions for the dissipative fluxes needs to be specified. In 
the present paper we assume that by the equilibration 
time Ti, the dissipative fluxes attained their longitudinal 
boost-invariant values. 
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(j}{x,p) is the deviation from equilibrium distribution 
function /'"■'. With shear viscosity as the only dissi- 
pative forces, (p{x,p) can be locally approximated by a 
quadratic function of 4-momcntum, 



(x,p) =£^,uP^p''. 



(4.3) 



Without any loss of generality e^^ can be written as 



as, 



^iiv = Ctt , C = — — ■ — - . 
2(e+p) 



(4.4) 



completely specifying the non-equilibrium distribution 
function. 



V. PARTICLE SPECTRA 

With the non-equilibrium distribution function thus 
specified, it can be used to calculate the particle spectra 
from the freeze-out surface. In the standard Cooper-Frye 
prescription, particle distribution is obtained as. 



E 



dN 



dN 



d^p dyd^pT 



d^^p^f{x,p) (5.1) 



In (r, x,y,rys) coordinate, the freeze-out surface is pa- 
rametcriscd as, 



E^ = {Tf{x,y)cosh'qs-,x,y,Tf{x,y)smh'qs), (5.2) 
and the normal vector on the hyper surface is, 



IV. NON-EQUILIBRIUM DISTRIBUTION 
FUNCTION 



With dissipation the system is not in equilibrium and 
the equilibrium distribution function. 



dTf dxf 
dT,^ ^ (cosh?7s, -- — , -- — , - smhr]s)Tfdxdydris 
dxf dyf 

. (^-^^ 
At the fluid position (T,x,y,rig) the particle 4- 

momenta are parameterised as. 



f^"Hx,p) 



exp[P{u^pt^ - ,u)] ± 1 ' 



(4.1) 



with inverse temperature (3 = l/T and chemical poten- 
tial ji can no longer describe the system. In a highly non- 
equilibrium system, distribution function is unknown. If 
the system is slightly off-equilibrium, then it is possible 
to calculate correction to equilibrium distribution func- 
tion due to (small) non-equilibrium effects. Slightly off- 
equilibrium distribution function can be approximated 



as. 



F{x,p)^f^"\x,p)[l + (b{x,p)l 



(4.2) 



p^ ~ {nnTCOsh{ris — Y),p^ ,p^ ,mTsinh(ris — Y)) (5.4) 
The volume element p^dY^f^ become, 

p^dT.fj_ — {mTCOsh{ri — Y)— px -V TTf)Tf dxdydrj (5.5) 

Equilibrium distribution function involve the term 
'-7^ which can be evaluated as, 



P^M/^ jimTCOshiji — y) — vt-Pt — fJ-fl) 
T ~ T 



(5.6) 



The non-equilibrium distribution function require the 



sum p^p^TT^^, 



PpPi^Tr'"' = aicosh^{ri — Y) + a2Cosh{ri — Y) + a^ (5-7) 

with 

ai = m2,(7r^^+T27r'"') (5.8) 

a2= -2mTip^n^'' + PyTT^y) (5.9) 

ag = p2^^:r ^ p2^TO + 2p^pyTT''y - m^j.r'^ tt'^\5 .10) 

I 



Inserting all the relevant formulas in Eg 15. II and inte- 
grating over spatial rapidity one obtains, 



dN 



dN'"i dN'""i 

+ 



dyd^PT dyd^pT dyd^px 



with. 



(5.11) 
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dN'""i 



dyd^PT (27r)3 



dN'"! 



dyd^PT (27r)^ 



dxdyTf[raTKi{nj3) — pT^TT/Ko^nf])] 



dxdyrfimTi^KsinP) + ^K2{np) + {^ + a-, + \)K,(np) + yKo(n/3)} 
-PT.VTT/{yK2(ri/3) + a2K^{n(3) + [^- + a-^ + l)i^o(n/3)}] 



r 



(5.12) 



(5.13) 



where Kq, Ki, K2 and K3 arc the modified Bessel func- 
tions. 

We will also show results for elliptic flow V2- It is de- 
fined as. 



Cl^^o<^^)d^ 



^2 - 27r rfjV 

Jo dyd'^pT 



(5.14) 



VI. RESULTS 
A. Evolution of the viscous fluid 

The energy-momentum conservation equations 
I19l:j.l()l:^.lll are solved using the SHASTA-FCT al- 
gorithm. We have made extensive changes to the 
publicly available code AZHYDRO (described in [3) 
for simulation of ideal fluid. The modifled code, called 
"AZHYDRO-KOLKATA", simulate the evolution of dis- 
sipative fluid, both in 1st and 2nd order theory. In this 
section we present results of AZHYDRO-KOLKATA for 
the 1st order theory of dissipative fluid. In the following 
we will show the results obtain in Au-|-Au collision at 
impact parameter b = 6.8 fm, which approximately 
corresponds to 16-24% centrality Au-t-Au collisions. 
With the same initial conditions, we have solved the 
energy-momentum conservation equations for ideal fluid 
and viscous fluid. 

In Fig.l, we have shown the constant energy density 
contour plot in x-y plane, after an evolution of 5 fm. The 
black lines are for ideal fluid evolution. The red and blue 
lines are for viscous fluid with ADS/CFT (77/3=0.08) and 
perturbative (77/3=0. 135) estimate of viscosity. Constant 
energy density contours, as depicted in Fig.l, indicate 
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FIG. 1: (color online). Contour plots of energy density at 
(proper) time r=5.6 fm. The black lines are for ideal fluid 
{■q/s—O). The red and blue lines are for viscous fluid with 
ADS/CFT (r7/s=0.08) and perturbative (77/3=0. 135) esti- 
mate of viscosity. 



that with viscosity fluid cools slowly. Cooling gets slower 
as viscosity increases. Thus at any point in the x-y plane, 
viscous fluid temperature is higher than that of the ideal 
fiuid. The results are in accordance with our expectation. 
For dissipative fluid, ideal equation of motion Ea l2.21l is 
changed to. 
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FIG. 2: (color online). Contour plots of x-component of fluid 
velocity v^ at r=5.6 fm. The black lines are for ideal fluid 
{rj/s=0). The red and blue lines are for viscous fluid with 
ADS/CFT and perturbative estimate of viscosity, 77/3=0. 08 
and 0.135. 



FIG. 3: (color online). Contour plots of temperature at y=0 
fm in a; — r plane. The black lines are for ideal fluid (r;/s=0). 
The red and blue lines are for viscous fluid with 77/5=0.08 and 
0.135 respectively. 



De^~{e + p)V^ 



7r^"'V<„Uj.> 



(6.1) 



Due to viscosity, evolution of energy density is slowed 
down. 

In Fig. 2, we have shown the constant v^ contour plot 
in x-y plane again at t=5.6 fm. As before the black lines 
are for the ideal fluid evolution. The red and blue lines 
are for viscous fluid with 77/5=0.08 and 0.135 respectively. 
In the central region of the fluid, viscous fluid has more 
velocity than its ideal counterpart. With viscosity while 
the energy density evolve slowly, the fluid velocity evolve 
faster. Contour plot of the y-component of fluid velocity 
also indicate similar results. 

To obtain an idea of transverse expansion of viscous 
fluid, as opposed to ideal fluid, in Fig. 3, we have shown 
the constant temperature contours in t — x plane , at a 
fixed value of y=0 fm. Transverse expansion is substan- 
tially enhanced in a viscous fluid. More the viscosity, 
more is the transverse expansion. The plot also indicate 
that at late time, fluid at x=y=0 behaves similarly to a 
ideal fluid. 

1st order dissipative theories are acausal. As men- 
tioned earlier, acausality can lead to unphysical behavior 
like reheating of the fluid in the early stage of evolu- 
tion |12. Ilq . Do we see any reheating? In FigQ] the 
evolution of temperature in viscous dynamics, with per- 
turbative estimate of viscosity (77/5=0.135) is shown. We 
have shown the temperature at two positions of the fluid. 



x=y=0 (the solid line) and x=0,y=3 fm (the dashed line). 
In both the positions of the fluid, with time as the fluid 
expands, temperature decreases (as it should be). We 
find no evidence of reheating. Reheating is not seen also 
with ADS/CFT estimate of viscosity. 

In Figl^l in 4-panels we have shown shear stress ten- 
sors TT'^'^{x,y = 0), T'^TT^'^{x,y = 0), 7r^^(a;,y = 0) and 
TTyy{x,y = 0) as a function of x. 77/5 =0.135. The solid, 
long-dashed, dashed and short-dashed lines are for time 
0.6, 2.2, 3.2 and 4.2 fm respectively. Initially at t=0.6 
fm, tt'^'^ is zero. As the fluid evolve, tt"^"^ increases rapidly 
to a maximum and then decreases. By 4 fm of evolution, 
it decreases to very small values. We also note that tt"^"^ 
is never very large. The viscous pressures t^tt'''', tt^^ 
and n^y are non-zero at initial time Ti=0.6 fm As the 
fluid evolve these viscous fluxes rapidly decreases to very 
small values. 

Viscosity generates entropy. In the model entropy gen- 
eration due to dissipation can be calculated as, 



d^S^' 



' fj,l^ 



277T 



(6.2) 



Evolution of spatially average entropy is shown in 
FigEl Entropy generation saturates after ^ 2 fm of evo- 
lution. It is expected. As seen in Fig|5l viscous fluxes 
rapidly decreases and by 2 fm of evolution, the viscous 
fluxes are decreased sufficiently and do not contribute 
significantly to the entropy. 
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FIG. 4; Evolution of the temperature in viscous dynamics 
with perturbative estimate of viscosity, ?7/s=0.135. The solid 
and dashed lines are for fluid at x=y=0 and x=0,y=3 fm 
respectively. The x=0,y=3 fm curve is plotted with the right 
side scale. 



FIG. 6; Evolution of average entropy with proper time, for 
two values of rj/s is shown. 



B. Particle spectra 
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FIG. 5: Shear stress tensor n'^'^{x,y = 0), T^n'^^{x,y = 0), 
n''''ix,y = 0) and 7r^«(j:,i/ = 0) at t=0.6, 2.2, 3.2 and 4.2 fm 
are shown in 4 panels. 



In this exploratory calculations wc have not attempted 
to fit experimental data. We just exhibit the effect of 
viscosity on (i) transverse momentum distribution and 
(ii) elliptic flow of pions. Viscosity influences the par- 
ticle production by (i) changing the freeze-out surface 
(freeze-out surface is extended) and (ii) by introducing a 
correction to the equilibrium distribution function. Non- 
equilibrium correction to equilibrium distribution func- 
tion depend, quadratically on the momentum and lin- 
early on the viscous fluxes. 

In Fig[7| we have shown the transverse momentum dis- 
tribution of pions obtained in the Cooper- Frye formalism. 
Freeze-out temperature is Tf=0.158 GeV. In this calcula- 
tion resonance contribution to pion spectra is neglected. 
Pion production is increased in viscous dynamics. We 
also note that effect of viscosity is more prominent at 
large px than at low px- Pt spectra of pions are flat- 
tened with viscosity. Particle production increases if vis- 
cosity increases. Thus while with ADS/CFT estimate of 
viscosity, 77/5=0.08, at pt =3 GeV, pion production is 
increased by a factor 3, with the perturbative estimate 
of viscosity, 77/5=0. 135, the production is increased by a 
factor of 5. Increase is even more at larger px- 

We have obtained the non-equilibrium distribution as 
a correction to the equilibrium distribution function. It 
is implied that non-equilibrium effects are small and the 
ratio 
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FIG. 7: Pt distribution of pions. The solid line is for ideal 
fluid. The long-dashed and medium-dashed lines are for 
viscous fluid with ADS/CFT (7;/s=0.Q8) and perturbative 
(77/s=0.135) estimate of viscosity. Non-equilibrium correction 
to equilibrium distribution function is included. 



FIG. 9: Elliptic flow as a function of transverse momen- 
tum. The solid line is for ideal fluid. The long-dashed, 
medium-dashed and short-dashed lines are for viscous fluid 
with 77/5=0. 08 and 0.135 respectively. Non-equilibrium cor- 
rection to equilibrium distribution function is included. 
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FIG. 8: Ratio of correction to particle production due to non- 
equilibrium distribution to equilibrium distribution function. 



is less than 1. In Fig|S| the ratio is shown as a function 
of PT- With ADS/CFT estimate of viscosity, rj/s^O.08, 
non-equihbrium correction to particle production become 
comparable to equilibrium contribution beyond Pt=^ 
GeV. However, with perturbative estimate, 77/5=6. 135, 
non-equihbrium correction become comparable to or ex- 



ceeds the equilibrium contribution at pT=4:.5 GeV. Thus 
with perturbative estimate of viscosity, hydrodynamic 
description break down above pr ^ i.bGeV. The blast 
wave model analysis 'l(l| on the otherhand indicated that 
viscous dynamics get invahdated beyond pT ^1.7 GeV. 
The results are not contradictory. In the blast wave 
model, at the freeze-out, viscosity is quite large, sound at- 
tenuation length Tg ~ 1.4 fm. In the present simulation, 
even for perturbative estimate of viscosity, sound atten- 
uation length at the freeze-out is Tg ~ 0.2fm, 7 times 
smaller than the sound attenuation length used in the 
blast wave analysis. Naturally, non-equilibrium correc- 
tions to equilibrium distribution function remains small 
over an extended px range. 

We have also calculated the elliptic flow in the model. 
Being a ratio, elliptic flow is very sensitive to the model. 
Experimentally, elliptic flow saturates at large px- It 
is known that ideal fluid does not explain the satura- 
tion of elliptic flow. In contrast to experiment, with 
ideal fluid, elliptic flow continues to increase with px- In 
Figl^l we have compared the elliptic flow in ideal and 
viscous fluid. The solid line is V2 for the ideal fluid. 
The long-dashed and medium-dashed lines are for vis- 
cous fluid with ADS/CFT (77/3=0.08) and perturbative 
(77/3=0.135) estimated viscosity. Elliptic flow decreases 
with viscosity. As viscosity increases, elliptic flow is also 
reduced. We also note that both for ADS/CFT and per- 
turbative estimate of viscosity, elliptic flow indicate sat- 
uration at large px- The result is very encouraging, as 
experimentally also elliptic flow tends to saturate at large 



10 



Pt- 

As discussed earlier, ideal fluid dynamics, can explain 
a large volume of data in Au+Au collisions at RHIC. Our 
present knowledge about the hot dense matter produced 
in central Au+Au collisions are obtained from the ideal 
fluid analysis. As shown in the present paper, QGP fluid, 
even with ADS/CFT estimate of viscosity ?7/s=0.08, gen- 
erate enough entropy to enhance particle production by 
a factor of 3 at PT—i GeV. Naturally, if QGP fluid is vis- 
cous, initial conditions as required to explain RHIC data 
with ideal fluid dynamics will overpredict the experimen- 
tal Pt distribution. Viscous fluid dynamics will require 
much less initial temperature than an ideal fluid to ex- 
plain the same pT spectra. As an example, in Fig. 10, we 
have compared the pion spectra obtained in viscous dy- 
namics with ADS/CFT estimate of viscosity (ry/s^O.OS), 
initialised with entropy density of 110,80, and 60 fm^^ 
with the pion spectra obtained in an ideal fluid dynamics, 
initialised with entropy density s=110 fm~^. For all the 
fluids, the initial time is ri=0.6 fm and freeze-out temper- 
ature is 158 MeV. Viscous fluid initialised with entropy 
density between 60-80 fm~^ , compare well with the pion 
spectra from ideal fluid initialised at much higher entropy 
density. To produce the same pion spectra, while ideal 
fluid require a initial temperature of 350 MeV, the vis- 
cous fluid require much less temperature between 270-290 
MeV. The ideal fluid dynamics can overestimate the ini- 
tial temperature of fluid produced in Au+Au collisions 
at RHIC by 20-30%. 
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FIG. 10: The solid circle is the pr distribution obtained in 
ideal fluid dynamics, with initial entropy density s=110/m~'^. 
The solid, long dashed and dashed lines are for viscous fluid 
with ADS/CFT estimate of viscosity, r;/s=0.08, initialised at 
entropy density s=60,80 and 110 fmT^ respectively. Non- 
equilibrium correction to equilibrium distribution function is 
included. 



VII. SUMMARY AND CONCLUSIONS 

We have studied the boost-invariant hydrodynamic 
evolution of QGP fluid with dissipation due to shear vis- 
cosity. In this study we have employed the 1st order 
theory of dissipative relativistic fluid. 1st order theories 
suffer from the problem of causality, signal can travel 
faster than light. Unphysical effects like reheating of 
the fluid, early in the evolution, can occur. However, 
for a fluid like QGP, where viscosity is small, with ap- 
propriate initial conditions, effects of causality violation 
can be minimised. In this model study, we have con- 
sidered two values of viscosity, the ADS/CFT motivated 
value, rj/s «0.08 and perturbatively estimated viscos- 
ity, rj/s «0.135. Both the ideal and viscous fluids are 
initialised similarly. At the initial time ri=0.6 fm, ini- 
tial central entropy density is 110 fm^^ , with transverse 
profile taken from a Glauber model calculation. Viscous 
hydrodynamics require initial conditions for the shear- 
stress tensor components. It is assumed that at the equi- 
libration time, the shear stress tensors components have 
reached their boost-invariant values. The initial condi- 
tions of the fluid are such that for both the values of 
viscosity (77/5=0. 08 and 0.135), the condition of valid- 
ity of viscous hydrodynamics, Fg/r << 1 is satisfied all 
through the evolution. Explicit simulation of ideal and 
viscous fluids confirms that energy density of a viscous 



fluid, evolve slowly than its ideal counterpart. The fluid 
velocities on the other hand evolve faster in viscous dy- 
namics than in ideal dynamics. Transverse expansion is 
also more in viscous dynamics. For a similar freeze-out 
condition freeze-out surface is extended in viscous fluid. 

We have also studied the effect of viscosity on par- 
ticle production. Viscosity generates entropy leading 
to enhanced particle production. Particle production 
is increased due to (i) extended freeze-out surface and 
(ii) non-equilibrium correction to equilibrium distribu- 
tion function. Non-equilibrium correction to equilibrium 
distribution function is a dominating factor influencing 
the particle production at large pr- With ADS/CFT 
(perturbative) estimate of viscosity, at pt=3 GeV, pion 
production is increased by a factor 3 (5) . Increase is 
even more at large px- While viscosity enhances particle 
production, it reduces the elliptic flow. At pr—^ GeV, 
for ADS/CFT(perturbative) estimate of viscosity, ellip- 
tic flow is reduced by a factor of 2(3). We also flnd that 
at large px elliptic flow tends to saturate. 

To conclude, present study shows viscosity, even if 
small, can be very important in analysis of RHIC Au+Au 
collisions. Currently accepted initial temperature of hot 
dense matter produced in RHIC Au+Au collisions, ob- 
tained from ideal fluid analysis can be changed by 20% 
or more with dissipative dynamics. 
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